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Abstract. Here results of numerical radiation hydrodynamical simula- 
' tions are presented which explore the energetic impact of massive stars 

on the interstellar medium. We study the evolution of the ambient 
^ . gas around isolated massive stars in the mass range between 15 and 

1 ^-^ ' 85 Mq in order to analyze the formation of structures and further the 

C/3 transfer and deposit of the stellar wind and radiation energy into the 

circumstellar medium until the stars explode as a supernovae. The 
derived energy transfer efficiencies are much smaller than analytically 
estimated and should be inserted into chemo-dynamical evolutionary 
models of galaxies as appropriate parameter values. As an additional 
issue the element release in the Wolf-Rayet phases and its detectability 
CO ■ have been investigated for comparison with observations. 

p 

1 Introduction 

o . 

. Besides external effects structure and evolution of galaxies is primarily determined 
by the energy budget of its interstellar medium (ISM). Reasonably, also the ISM 
structure and both its radiative and kinetic energy contents are determined by the 
k> " energy deposit of different sources, on the one hand, and by the energy loss by 
. means of radiative cooling on the other. Mainly massive stars contribute signifi- 
ed ' cantly to this structure formation like e.g. cavities, holes, and chimneys in the Hi 
gas and superbubbles of hot gas (Recchi & Hensler 2006). On large scales the en- 
ergy release by massive stars triggers the matter circulations via galactic outflows 
from a gaseous disk and galactic winds (Recchi & Hensler, this volume). By this, 
also the chemical evolution is affected thru the loss of metal-enriched gas from a 
galaxy (see e.g. Recchi et al. I2006p . 

Supernova (SN) explosions as an immediate consequence of star formation (SF) 
stir up the ISM by the expansion of hot bubbles, deposit turbulent energy into 
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the ISM and, thereby, regulate the SF again (Hensler & Rieschick [2002) . This 
negative energy feedback is enhanced at low gravitation because the SN energy 
exceeds easily the galactic binding energy and drives a galactic wind. Vice versa, 
SN and stellar wind-driven bubbles sweep up surrounding gas and can, by this, 
excite SF self-propagation as a positive feedback mechanism (e.g. Ehlerova et al. 
fT997l Fukuda & HanawaUDDO]). 

Since the numerical treatment of the galactic chemo- dynamical evolution, al- 
though when it deals with two gas phases, cannot resolve the ISM spatially suffi- 
ciently, for a detailed consideration of small-scale processes the chemo-dynamical 
modeling has to apply parametrizations of plasmaphysical processes. It should, 
however, be repeatedly emphasized that an appropriate chemo-dynamical de- 
scription cannot vary those parameters arbitrarily, because they rely on results 
from theoretical, numerical, and/or empirical studies and are by this strongly 
constrained. This has been described and applied to different levels of chemo- 
dynamical models (for an overview see Hensler, 120031 and references therein, and 
more recently Hensler et al. 120041 Harfst et al.. l2006p . 

As an example, analytical consideration of heating and cooling balance by 
Koppen et al. (|1995p making use of the cooling by coUisionally excited radiation 
(Bohringer & Hensler I1989P yields that the SF rate under reasonable ISM states 
is dependent on the square of the gas density with a slight influence of the heating 
efficiency. This fact is surprising: although the ionizing radiative fluxes of massive 
stars are pretty well known from static stellar atmosphere models (see e.g. Thomp- 
son [1984]) or more sophisticated wind atmospheres (Kudritzki & Puis 12000]) . the 
fraction of the radiative energy that is thermalizing the gas is uncertain. 

In addition, massive stars act also dynamically on their surrounding ISM by 
a strong stellar wind (Kudritzki & Puis 12000]) . Although the stellar wind power 

= ^MyjV^ with the mass-loss rate (Mw) and the terminal wind velocity {vw) 
can be easily evaluated from model and observational values, the fraction of the 
wind luminosity which is transferred e.g. into thermal energy, the so-called thermal 
energy transfer efficiency, or into turbulent energy is not obvious from principles. 

2 Theory of radiative and stellar wind bubbles 

Within the framework of the simplest theoretical approach an O star suddenly 
"turns on" in a constant density medium at rest and begins to ionize its surround- 
ings. The reaction of the medium to the stellar photons is well-known and has 
been described in detail in standard textbooks. With the particular knowledge of 
the time evolution of the ionization front at distance from a star of age r and with 
respect to the initial Stromgren radius Rs which depends on the Lyman contin- 
uum photon luminosity, the Hi number density of the ambient medium, and the 
isothermal sound speed of the ionized gas Lasker ([1967P derived the kinetic energy 
of the expanding swept-up shell E]^, the ionization energy stored in the Hii region 
Ei, and the thermal energy of warm gas in the Hii region Et (see also Freyer et 
al. I2003|, hereafter: Paper I). As an appropriate approximation the temperature 
in the Hii region needed in the set of formulae can be set to 8000 K. 
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Also for stellar wind bubbles (SWBs) the radial structure can be analytically 
derived (see Weaver et al. 119771 and Paper I) and, furthermore, also the kinetic 
energy as E]^ = L^t and the thermal energy of the hot gas as i^t = ^ Lw t. 
Besides heat conduction and cooling in the hot bubble here also neglected is the 
fact that part of the thermal energy might be used for coUisional ionization of gas. 
Thus, the results for and Et are upper limits. 

Because these analytical derivations cannot account for aspherical structure 
formation and turbulence within the SWB which both consume stellar energy 
and alter the energy fractions numerical models are needed to investigate the 
massive stellar energy deposit into the ISM. For a comprehensive tabulation and 
comparison of former models the interested reader is referred to Paper I. 

3 The model 

The numerical code applied for this purpose and used to obtain the results pre- 
sented here is described in Paper I. The hydrodynamical equations are solved 
together with the transfer of H-ionizing photons on a 2d cylindrical grid. The time- 
dependent ionization and recombination of hydrogen is calculated each timestep, 
and we carefully take stock of all the important energy exchange processes in the 
system. A comprehensive description of the algorithm is given in Yorke & Kaisig 
(fTMS]) . Yorke & Welz (fT996l) . and in Paper I. 

This code is applied to the evolution of massive stars with 4 different masses: 
15 Mq (Kroeger, Freyer, Hensler & Yorke, in prep, as Paper IV of this series), 
35 Mq (Freyer et al. ,200 6: Paper II of this series), 60 Mq (Paper I), and 85 Mq 
(Kroeger et al. I2006bl Paper III of this series). Because we started the project 
also at first for comparison with the 60 Mq model published by Garcia-Segura 
et al. |l996a) we adopted basically the same stellar parameters applied by them 
from stellar evolutionary models and from observations. As initial conditions of 
the ambient medium we also use the same undisturbed background gas layer as 
Garcia-Segura et al. with Hi number density ng = 20 cm~^ and temperature Tq — 
200 K. This yields a thermal pressure "typical" for the ISM at the solar galacto- 
centric distance. Also the 35 Mq model was aimed for comparison with Garcia- 
Segura et al. J1996b|) and used their stellar parameters along its life. Because 
from this group no evolutionary tracks and stellar parameters of other stellar 
masses became available yet, however, for the 15 Mq and the 85 Mq models we 
were compelled to apply stellar evolutionary models from Schaller et al. (I1992p . 
Although this makes the sample inhomogeneous, the extensive computational time 
necessary for the models only allowed until now comparison of the two 60 Mq star 
model sources. As a first result, we obtained for the Schaller et al. model a by one 
order of magnitude weaker wind luminosity (Hensler, Kroeger, & Freyer, in prep.). 
Furthermore, the evolutionary sequence of Wolf-Rayet (WR) and Luminous Blue 
Variable (LBV) phases are reverse from one group to the other in the sense that 
the WR phase precedes the LBV phase in the Langer models. 

In our models the stellar parameters Mu,, v^, effective temperature (Tcff), and 
photon luminosity in the Lyman continuum (iLyc) are time-dependent boundary 
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Fig. 1. Density plot of the wind-blow and radiation-driven bubble around a 60 Mq 
main-sequence star at an age of 140 000 years. 



conditions which drive and govern the evolution of the circumstellar gas. 

We start our calculations with the sudden turn-on of the zero-age main-sequence 
stellar radiation field and stellar wind. For the surrounding gas we neglect any 
photo-dissociation of molecular content here which would in reality precede the 
photo-ionization stage (e.g. Kroeger, Hensler, & Freyer, in prep, as Paper V of 
this series). 

4 Results 

As a first result, the winds structure the surrounding gas by dynamical shock 
front instabilities during its vehement expansion. These lead to the formation 
of condensations within the wind shell which, on the other hand, have higher 
optical depth for the ionizing photons. Finger- like extentions of the Hii region 
surrounding the SWB are most characteristicly formed for the 60 M© model (see 
Fig.l) from photons penetrating deeper into the ambient ISM. In addition, photo- 
evaporation peels off the inner boundary of the wind shell and drives mass loading 
and turbulence. 

Another relevant question addressedto these models is their energy deposit 
into the ISM. Here we can only briefly mention a few major issues. Over the 
whole evolution one can trace the energy released in different forms as kinetic i?k, 
photo- ionizing Ei, and thermal Et energies, respectively, by book keeping so that 
at every stellar age r each integrated energy form released can be compared with 
that inherently available in the computational domain. The ratio of deposited 
stellar energy to the increase in one form compared to the initial state can be 
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Fig. 2. Energy transfer efficiencies e^, £i, and et for kinetic, ionizing photon, and thermal 
energies for the 15, 35, 60, and 85 M0 models. k,windless and k,wind as well as for i and 
t, respectively, represent the analytical values according to egu M.ll and its extention by 
Lw, while ek, ei, et give the efficiency values at the end of the stellar lifetimes. 



defined as an energy transfer efficiency e. As the simplest approach, e for pure 
Lyman continuum radiation into the other energy forms is then descibed by the 
analytical approach of 

where E can be any of i?k, ^'i, and Et, respectively. For radiation plus wind energy 
this nominator has to be extended accordingly to ( (£Lyc) + (Lw))- 

These e values have been derived for all four massive star models and their 
energy forms (Hensler, Freyer, & Kroeger, 2007, in prep, as Paper VI of this 
series) and are exhibited in Fig. 2. Not surprisingly, the windless values are below 
the wind models except the ionization efficiency ei. Nevertheless, the average 
e only exceed one percent for the thermal efficiency of the 60 Mq wind model. 
As already mentioned, there the structure formation is strongest and might be 
reduced applying the Schaller et al. (|1992p model. 

The unexpected decrease of ei (the upper dashed curve in Fig. 2) with increasing 
stellar mass is caused by the growing strength of stellar winds with mass that have 
the dynamical effect to compress the ambient shells much stronger so that the 
recombination is enhanced. This means, ionization energy is more efficiently lost 
for larger stellar masses. 

Comparing the windless i.e. pure photo-ionizing models with analytical results 
e.g. by Lasker (|1967p of about one percent, our efficiencies fall short by more than 
one order of magnitude even to below 0.1 percent (Fig. 2). et reaches more than 



118 



CRAL-2006. Chemodynamics: From First Stars to Local Galaxies 



one percent only for 60 M0 . As another important issue, it should be mentioned 
that the final efficiencies {ek,ei,et in Fig. 2) are located significantly below the 
values integrated over the stellar life because of the waeker wind at the last stages. 

The efficiencies derived from our models are much less than expected and of 
fundamental importance for the energy budget of galaxies during their evolution. 
They should be easily inserted into chemo-dynamical evolutionary codes because 
integrating the stellar energy release by massive stars over their concedingly short 
lifetimes without taking detailed heating vs. cooling processes into account would 
falsify the treatment of the ISM significantly. A comprehensive discussion is in 
preparation as Paper VI. 

Finally, a byproduct of these studies should be mentioned, namely, that the 
elements release in the WR phase of most massive stars, in particular of carbon, 
although at first incorporated into the hot SWB gas are enabled to cool by mixing 
effects and become observable as warm ionized Hii gas (Kroeger et al. I2006ap . 

The author thanks Danica Kroeger and Tim Freyer for their strong engagement to this project 
and the organizers for this gorgeous and most exciting conference. 
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